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O . ABSTRACT 

(N 

We assume that the early universe is homogeneous, anisotropic, 
Q . and is dominated by the mutually BPS 22'55' intersecting branes 

O ! of M theory. The spatial directions are all taken to be toroidal. 

Using analytical and numerical methods, we study the evolution 
of such an universe. We find that, asymptotically, three spatial 
^ ■ directions expand to infinity and the remaining spatial directions 

j^ . reach stabilised values. Any stabilised values can be obtained by 

a fine tuning of initial brane densities. We give a physical descrip- 
tion of the stabilisation mechanism. Also, from the perspective of 
four dimensional spacetime, the effective four dimensional New- 
ton's constant G4 is now time varying. Its time dependence will 
follow from explicit solutions. We find in the present case that, 
asymptotically, G4 exhibits characteristic log periodic oscilla- 
tions. 



PACS numbers: 11.25.Yb, 98.80.Cq, ll.25.-w 



I. Introduction 



In early universe, temperatures and densities reach Planckian scales. Its 
description then requires a quantum theory of gravity. A promising candidate 
for such a theory is string/M theory. When the temperatures and densities 
reach string/M theory scales, the appropriate description is expected to be 
given in terms of highly energetic and highly interacting string/M theory 
excitations [T] - [S]. u In this context, one of us have proposed in an earlier 
work an entropic principle according to which the final spacetime configura- 
tion that emerges from such high temperature string/M theory phase is the 
one that has maximum entropy for a given energy. This principle implies, 
under certain assumptions, that the number of large spacetime dimensions 

is 3 + 1 [g. 

High densities and high temperatures also arise near black hole singular- 
ities. Therefore, it is reasonable to expect that the string/M theory config- 
urations which describe such regions of black holes will describe the early 
universe also. 

Consider the case of black holes. Various properties of a class of black 
holes have been successfully described using mutually BPS intersecting con- 
figurations of string/M theory branes. u Black hole entropies are calculated 
from counting excitations of such configurations, and Hawking radiation is 
calculated from interactions between them. 

In the extremal limit, such brane configurations consist of only branes 
and no antibranes. In the near extremal limit, they consist of a small num- 
ber of antibranes also. It is the interaction between branes and antibranes 
which give rise to Hawking radiation. String theory calculations are tractable 
and match those of Bekenstein and Hawking in the extremal and near ex- 
tremal limits. But they are not tractable in the far extremal limit where the 
numbers of branes and antibranes are comparable. However, even in the far 
extremal limit, black hole dynamics is expected to be described by mutually 



^ Only string theory is considered in these references. But their arguments can be 
extended for M theory also leading to similar conclusions. 

^ Mutually BPS intersecting configurations means, for example, that in M theory two 
stacks of 2 branes intersect at a point; two stacks of 5 branes intersect along three common 
spatial directions; a stack of 2 branes intersect a stack of 5 branes along one common spatial 
direction; waves, if present, will be along a common intersection direction; and each stack 
of branes is smeared uniformly along the other brane directions. See [9] for more details 
and for other such string/M theory configurations. 



BPS intersecting brane configurations where they now consist of branes, an- 
tibranes, and other excitations hving on them, all at non zero temperature 
and in dynamical equilibrium with each other [10] - [19]. For the sake of 
brevity, we will refer to such far extremal configurations also as brane con- 
figurations even though they may now consist of branes and antibranes, left 
moving and right moving waves, and other excitations. 

The entropy S" of N stacks of mutually BPS intersecting brane configura- 
tions, in the limit where S" S> 1 , is expected to be given by 
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5* ~ ]^ y/nj + nj ~ £ 2 (1) 

where nj and nj , / = 1, ■ ■ ■ , A^ , denote the numbers of branes and an- 
tibranes of /*'' type, S is the total energy, and the second expression applies 
for the charge neutral case where n/ = nj for all / . The proof for this 
expression is given by comparing it in various limits with the entropy of the 
corresponding black holes [IOl[TT], see also [12] - [21]. For A^ < 4 and when 
other calculable factors omitted here are restored, this expression matches 
that for the corresponding black holes in the extremal and near extremal 
limit and, in the models based on that of Danielsson et al [12], matches upto 
a numerical factor in the far extremal limit [10] - [21] also. However, no such 
proof exists for A^ > 4 since no analogous object, black hole or otherwise, is 
known whose entropy is oc S* with * > 2 . 

Note that, in the limit of large S , the entropy S{S) is '^ S for radiation in 
a finite volume and is ^ S for strings in the Hagedorn regime. In comparison, 
the entropy given in ([1]) is much larger when N > 2 . This is because the 
branes in the mutually BPS intersecting brane configurations form bound 
states, become fractional, and support very low energy excitations which 
lead to a large entropy. Thus, for a given energy, such brane configurations 
are highly entropic. 

Another novel consequence of fractional branes is the following. Accord- 
ing to the 'fuzz ball' picture for black holes [22], the fractional branes arising 
from the bound states formed by intersecting brane configurations have non 
trivial transverse spatial extensions due to quantum dynamics. The size of 
their transverse extent is of the order of Schwarzschild radius of the black 
holes. Therefore, essentially, the region inside the 'horizon' of the black hole 
is not empty but is filled with fuzz ball whose fuzz arise from the quantum 
dynamics of fractional strings/branes. 

Chowdhury and Mathur have recently extended the fuzz ball picture to 



the early universe [201 [2T]. They have argued that the early universe is 
filled with fractional branes arising from the bound states of the intersecting 
brane configurations, and that the brane configurations with highest A^ are 
entropically favorable, see equation ([1]). 

However, as mentioned below equation ([1]) and noted also in [2D1|2T], the 
entropy expression in ([T]) is proved in various limits for A^ < 4 only and 
no proof exists for A^ > 4. Also, we are not certain of the existence of any 
system whose entropy S{S) is parametrically larger than S^ for large S . 
See related discussions in [231 121] • Therefore, in the following we will assume 
that A^ < 4 . Then, a homogeneous early universe in string/M theory may be 
taken to be dominated by the maximum entropic A^ = 4 brane configurations 
distributed uniformly in the common transverse space. 

Such A^ = 4 mutually BPS intersecting brane configurations in the early 
universe may then provide a concrete realisation of the entropic principle 
proposed earlier by one of us to determine the number (3 + 1) of large space- 
time dimensions [8]. Indeed, in further works [231 1211 125], using M theory 
symmetries and certain natural assumptions, we have shown that these con- 
figurations lead to three spatial directions expanding and the remaining seven 
spatial directions stabilising to constant sizes. 

In this paper, we assume that the early universe in M theory is homo- 
geneous and anisotropic and that it is dominated by A^ = 4 mutually BPS 
intersecting brane configurations, g In this context, it is natural to assume 
that all spatial directions are on equal footing to begin with. Therefore we as- 
sume that the ten dimensional space is toroidal. We then present a thorough 
analysis of the evolution of such an universe. 

The corresponding energy momentum tensor Tab has been calculated 
in [20] under certain assumptions. However, general relations among the 
components of Tab niay be obtained [24j using U duality symmetries of M 
theory which are, therefore, valid more generally, u We show in this paper 
that these U duality relations alone imply, under a technical assumption, that 
the A^ = 4 mutually BPS intersecting brane configurations with identical 
numbers of branes and antibranes will asymptotically lead to an effective 
(3 + 1) dimensional expanding universe. 

In order to proceed further, and to obtain the details of the evolution. 



^ There is an enormous amount of work on the study of early universe in string/M 
theory. For a smah, non exhaustive, sample of such works, see [26] - [38] , 

* Such U duality relations are present in the case of black holes also. We point them 
out in Appendix A. 



we make further assumptions about T^b ■ We then analyse the evolution 
equations in D dimensions in general, and then specialise to the eleven di- 
mensional case of interest here. 

We are unable to solve exphcitly the relevant equations. However, apply- 
ing the general analysis mentioned above, we describe the qualitative features 
of the evolution of the A^ = 4 brane configuration. In the asymptotic limit, 
three spatial directions expand as in the standard FRW universe and the 
remaining seven spatial directions reach constant, stabilised values. These 
values depend on the initial conditions and can be obtained numerically. 
Also, we find that any stabilised values may be obtained, but requires a fine 
tuning of the initial brane densities. 

Using the analysis given here, we present a physical description of the 
mechanism of stabilisation of the seven brane directions. The stabilisation 
is due, in essence, to the relations among the components of Tab which 
follow from U duality symmetries, and to each of the brane directions in the 
A^ = 4 configuration being wrapped by, and being transverse to, just the 
right number and kind of branes. This mechanism is very different from the 
ones proposed in string theory or in brane gas models [2B] - [21] to obtain 
large 3 + 1 dimensional spacetime. (See section I A below also.) 

In the asymptotic limit, the eleven dimensional universe being studied 
here can also be considered from the perspective of four dimensional space- 
time. One then obtains an effective four dimensional Newton's constant G4 
which is now time varying. Its precise time dependence will follow from 
explicit solutions of the eleven dimensional evolution equations. 

We find that, in the case of A^ = 4 brane configuration, G4 has a char- 
acteristic asymptotic time dependence : the fractional deviation 6G4 of G4 
from its asymptotic value exhibits log periodic oscillations given by 

SG4 oc — Sin{u Int + cf)) . (2) 

The proportionality constant and the phase angle ^ depend on initial condi- 
tions and matching details of the asymptotics, but the exponents a and u 
depend only on the configuration parameters. Such log periodic oscillations 
seem to be ubiquitous and occur in a variety of physical systems [391 SQlIlT]. 
But, to our knowledge, this is the first time it appears in a cosmological 
context. One expects such a behaviour to leave some novel imprint in the 
late time universe, but its implications are not clear to us. 



Since we are unable to solve the evolution equations explicitly, we anal- 
yse them using numerical methods. We present the results of the numerical 
analysis of the evolution. We illustrate the typical evolution of the scale fac- 
tors showing stabilisation and the log periodic oscillations mentioned above. 
By way of illustration, we choose a few sets of initial values and present the 
resulting values for the sizes of the stabilised directions and ratios of the 
string/M theory scales to the effective four dimensional scale. 

We also discuss critically the implications of our assumptions. As we will 
explain, many important dynamical questions must be answered before one 
understands how our known 3-1-1 dimensional universe may emerge from 
M theory. Until these questions are answered and our assumptions justified, 
our assumptions are to be regarded conservatively as amounting to a choice 
of initial conditions which are fine tuned and may not arise naturally. 

The organisation of this paper is as follows. In section II, we describe 
the U duality symmetries of M theory and their consequences, and present 
our ansatzes for the energy momentum tenosr Tab and for the equations 
of state. In section III, we present a general analysis of D dimensional 
evolution equations. In section IV, we specialise to the eleven dimensional 
case of A^ = 4 intersecting brane configurations and describe the various 
results mentioned above. In section V, we discuss the stabilised values of the 
brane directions, their ranges, and the necessity of fine tuning. In section 
VI, we present the four dimensional perspective and the time variations 
of G4 . In section VII, we present the results of numerical analysis. In 
section VIII, we conclude by presenting a brief summary, a few comments 
on the assumptions made, and by mentioning a few issues which may be 
studied further. In Appendix A, we highlight the points related to U duality 
symmetries in the black hole case. In Appendices B - D, we present certain 
results required in the text of the paper. 

A. Intersecting brane vs Brane gas models 

In this subsection, we note that the branes and antibranes in the mutu- 
ally BPS intersecting brane configurations considered here and in [20] - [25] 
are different from those in the string/brane gas models [28] - [31] in many 
important aspects. These differences are explained in detail in section 2.6 of 
[20] and section 6 of [21] ■ Briefly, the differences are the following. 

(1) In brane gas models, the branes can intersect each other arbitrarily. In the 
brane configurations here, the intersections must follow specific rules. Con- 



sequently, U duality symmetries of M theory imply certain relations among 
the components of the energy momentum tensor Tab which turn out to be 
crucial elements in our case |23l [2] 



(2) The branes in brane gas models support excitations on their surfaces and, 
at high energies, have S ^ S where S is the entropy and £ the energy. Here, 
the intersecting branes form bound states, become fractional, support very 
low energy excitations and, hence, are highly entropic. At high energies, S ~ 

N 

S~ which, for N > 2 , vastly exceeds the entropy in brane gas models. Such 
intersecting brane configurations are, therefore, the entropically favourable 
ones. 

(3) In brane gas models, the branes are assumed to annihilate if they inter- 
sect each other. Here, the intersections are necessary for formation of bound 
states, and thereby of fractional branes leading to high entropic excitations. 
Also, the intersecting configurations here consist of branes, antibranes, and 
large number of low energy excitations living on them. All these constituents 
are at non zero temperature and in dynamical equilibrium with each other 
[in]-[19]. The excitations are long-lived and non interacting to the lead- 
ing order, hence the branes and antibranes here are metastable and do not 
immediately annihilate. 

Note that, in string/M theory description, the Hawking evaporation of 
black holes is due to the annihilation of branes and antibranes. Also, large 
black holes consist of stacks of large numbers of branes and antibranes, and 
have a long life time. This implies that such stacks of branes and antibranes 
do not immediately annihilate and have a long life time. Hence, we assume 
that the mutually BPS intersecting brane configurations describing the early 
universe also consists of stacks of large numbers of branes and antibranes 
having a long life time and, in particular, that the brane antibrane annihi- 
lation effects are negligible during the evolution of the universe atleast until 
the brane directions are stabilised resulting in an effective 3 + 1 dimensional 
universe. 



II. U duality symmetries and Equations of state 



In this paper, we assume that the early universe in M theory is homo- 
geneous and anisotropic and that it is dominated by A^ = 4 mutually BPS 



intersecting brane configurations. To be specific, we consider 22'55' configu- 
rations, u We study the consequent evolution of sucli an universe dictated 
by a (10 + 1) dimensional effective action given, in tlie standard notation, by 



Sn = j^r^:r- d''x./^R + Str (3) 

iOTrCrii J 

wliere S^r is tlie action for tlie fields corresponding to the branes. The cor- 
responding equations of motion are given, in the standard notation and in 
units where SttGh = 1 , by |j 



Rab - I gABR = Tab , ^ V^T^^ = (4) 

A 
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where A = (0, i) with i = 1, 2, ■ ■ ■ , 10 and T^b is the energy momentum 
tensor corresponding to the action 5';,^ . 

For black hole case. Tab is obtained from the action for higher form 
gauge fields. With a suitable ansatz for the metric, equations of motion 
(HI) are solved to obtain black hole solutions. For cosmological case. Tab is 
often determined using equations of state of the dominant constituent of the 
universe. Such equations of state may be obtained if the underlying physics 
is known; or, one may assume a general ansatz for them and proceed. LJ 

Tab for intersecting branes in the early universe has been calculated in 
[20] assuming that the branes and antibranes in the intersecting brane config- 
urations are non interacting and that their numbers are all equal, i.e. rij = nj 
for / = 1, 2, ■ ■ ■ , A^ and ni = ■ ■ ■ = un . However, general relations among 
the components of Tab inay be obtained |23] using U duality symmetries of 
M theory, involving chains of dimensional reduction and uplifting and T and 
S dualities of string theory, using which 2 branes and 5 branes or 22'55' and 



^ In our notation, 22'55' denotes two stacks each of 2 branes and 5 branes, all inter- 
secting each other in a mutually BPS configuration. Similarly for other configurations, 
e.g. 55'5"W denotes three stacks of 5 branes intersecting in a mutually BPS configuration 
with a wave along the common intersection direction. 

^ In the following, the convention of summing over repeated indices is not always 
applicable. Hence, we will always write the summation indices explicitly. Unless indicated 
otherwise, the indices A,B,--- run from to 10, the indices i,j, ■ ■ ■ from 1 to 10, and the 
indices I,J,--- from 1 to TV. 

^ This is similar to the FRW case. Equation of state P = f for radiation, or p = for 
pressureless dust, may be obtained from the physics of radiation or of massive particles; 
or, one may assume a general ansatz p = wp and proceed. 



55'5"W configurations can be interchanged. Such relations are vahd more 
generally, for example even when rij and nj are all different. 

These general relations on the equations of state are sufficient to show, 
under a technical assumption, that the A^ = 4 mutually BPS intersecting 
brane configurations with identical numbers of branes and antibranes, i.e. 
with 77-1 = ■■■ = 714 and ni = ■ ■ ■ = n^ , will asjTiiptotically lead to an 
effective (3 + 1) dimensional expanding universe. To obtain the details of the 
evolution, however, we need further assumptions and an ansatz of the type 
p = wp pit [25]. 

We now present the details. Let the line element ds be given by 

ds^ = -dt^ + Y^e^>^\dx'y (5) 

i 

where e^' are scale factors and, due to homogeneity. A* are functions of the 
physical time t only. (Parametrising the scale factors as e^* turns out to be 
convenient for our purposes.) It then follows that Tab depends on t only and 
that it is of the form 

T^B = diag{-p, p,) . (6) 

We assume that p > . From equations (jl]) one now obtains 

A?-E(Ai)' = 2p (7) 

i 

Xl, + AtXl = P^ + l{p-T.P^) (8) 

Pt + pKt + Y.P^K = (9) 

i 

where A = J^i^^ and the subscripts t denote time derivatives. Note, from 
equation ([7]), that Aj cannot vanish. Hence, if Aj > at an initial time to 
then it follows that e^ increases monotonically for t > to . We assume the 
evolution to be such that e"^ — )• oo eventually. 

In the context of early universe in M theory, it is natural to assume that 
all spatial directions are on equal footing to begin with. Therefore we assume 
that the ten dimensional space is toroidal. Further, we assume that the early 
universe is homogeneous and is dominated by the 22'55' configuration where, 
with no loss of generality, we take two stacks each of 2 branes and 5 branes 
to be along {x^,x^) , {x^,x^) , {x^,x^,x^,x^,x'^) , and {x"^ , x"^ , x^ , x^ , x"^) 
directions respectively, and take these intersecting branes to be distributed 



uniformly along the common transverse space directions {x^,x^,x^^) . Note 
that the total brane charges must vanish, i.e. rij = fij for all / , since the 
common transverse space is compact. We denote this 22'55' configuration 
as (12,34, 13567,24567) . The meaning of this notation is clear and, below, 
such a notation will be used to denote other configurations also. 

A. U duality relations 

We now describe the relations which follow from U duality symmetries, 
involving chains of dimensional reduction and uplifting and T and S dualities 
of string theory. See [2^j for more details. Let J,^ and tfc denote dimensional 
reduction and uplifting along k^^ direction between M theory and type IIA 
string theory, Ti denote T duality along z*^ direction in type IIA and IIB 
string theories, and S denote S duality in type IIB string theory. Then 
U dualities of the type tj TiSTi J,^ interchange i and j directions, and U 
dualities of the type '1'^ TjTj \,k transform one mutually EPS N intersecting 
brane configuration to another. 

For example, the U duality ts T^T^ is transforms a 2 brane configura- 
tion (12) to the 5 brane configuration (12345) . Similarly, the U duality 
t5 T1T2 Is transforms the 22'55' configuration (12, 34, 13567, 24567) to the 
Wbb'b" configuration (5, 12345, 23567, 14567) ; whereas te ^47^5 U trans- 
forms it to the 5'2'52 configuration (12456, 35, 13467, 27) . 

The U dualities transform the corresponding gravitational fields also. In 
the case of metric of the form given in equation ([5]), the U duality tfc TiTj J,^ 
transforms the A*s in the scale factors to A'*s given by 

A'^ = A^' - 2A , A'^' = A* - 2A , A'^ = A^ - 2A 

A'' = A' + A, /^z,j,fc, A= ^^ + ^J + ^' . (10) 

Consider the 2 brane configuration (12) with scale factors e^'' and the 
5 brane configuration (12345) with scale factors e'^ ' . By inspection, or by 
using U dualities tj T^STi \,j for appropriate (i, j), these scale factors may 
be expected to obey the 'obvious' symmetries: 

2 : A^ = A^ , A3 = . . . = A^o (11) 

5 : A'i = --- = A'^ , A'6 = --- = A'i° . (12) 

Now, these two configurations are related by the U duality ts ^^324 is . 
Hence the corresponding A*s and A'*s must obey the relations of the type 

10 



given in fllOp . Combined with the obvious symmetry relations above, it is 
straightforward to show that 

A" + 2A^ = 2A'II + A'^ = (13) 

where the superscripts || and _L denote spatial directions along and trans- 
verse to the branes respectively. 

Similarly, the obvious symmetry relations for the 22'55' configuration 
(12,34,13567,24567) are 

22'55' : A^ = A'^ = A' , A^ = A^ = A^^ . (14) 

Proceeding as in the case of 2 and 5 branes above, and using the U duality 
'1^5 T1T2 4,5 which relates the 22'55' and VF55'5" configurations, one obtains 
two more relations given by 



Ai + A^ + A^ = A2 + A=^ + A^ = . (15) 

In general, for an A^ intersecting brane configuration, the U duality sym- 
metries will lead to 10 — iV relations among the A*s, leaving only A^ of them 
independent. These relations are of the form X^i qA* = where q are con- 
stants. Clearly, such a relation can be violated by constant scaling of x* 
coordinates. Hence, we interpret it as implying a relation among the com- 
ponents of Tab ■ In view of equation ([S]), we interpret a U duality relation 
Y,i CjA* = as implying that 

^Qf = , r = P^ + \{p-Y.P3) ■ (16) 

i 3 

Substituting equation (ITB]) in equation (jH]), it follows upon an integration 

that 

^QAi = ce-^ (17) 

i 

where c is an integration constant. If Yl,iCi\\ = initially at t = to then 
c = . In such cases then J^i^iK = ^^ ^^^ ^ ^^d, hence, J^i^iX'' = v 
where v is another integration constant. 

In general J2i Ci\\ 7^ initially at t = to and, hence, c 7^ . Let the 
evolution be such that e^ ~ t^ — > 00 in the limit t — > 00 . Then it follows 
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from equation fll7p that J2i CiK — ^ in this hmit. If, furthermore, /3 > 1 
then equation ( 1T7|) also gives 

Y^CiX' = v + 0{t^-^)^v (18) 



where v is an integration constant. If /? < 1 then J2i c^A* is a function of t . 
We will see later that, in the solutions we obtain with further assumptions, 
P turns out to be > 1 for A^ > 1 . 

Note that, as can be seen from the above steps, the U duality relations 
follow as long as the directions involved in the U duality operations are 
isometry directions. Since none of the common transverse directions are 
involved in obtaining the relations above, it follows that they are valid even 
if the common transverse directions are not compact. Thus the U duality 
relations are applicable in such cases also. 

Similarly, the time dependence of A*s played no role in obtaining the U 
duality relations here. Hence, these relations may be expected to arise for 
the black hole case also. They indeed arise as we point out in Appendix A. 

B. A general result 

We now consider a general result for the 22'55' configuration that follows 
from the U duality relations alone [25j. The A*s for this configuration obey 
the relations given in equations flMl) and flTSl) . Note that a suitable U duahty, 
for example te T^T^ Iq , can transform 2 branes and 5 branes into each other. 
Hence, we will refer to two types of branes as being identical if they have 
identical numbers of branes and antibranes, i. e. P^ type is identical to J*'* 
type if ni = nj and nj = nj . 

Consider the case when 2 and 2' branes in the 22'55' configurations are 
identical. This will enhance the obvious symmetry relations. It is easy to 
see that we now have one more independent relation A^ = A'^ . If, instead, 5 
and 5' branes are identical, then the extra independent relation is A^ = A^ . 
Similarly, if 2 and 5' branes are identical then, after a few steps involving U 
duality te T4T5 Iq which interchanges 2 and 5' branes, it follows that the 
extra independent relation is A^ = A^ . 

Now if all the four types of branes in the 22'55' configuration are identical, 
i.e. if ni = ■ ■ ■ = 77.4 and ni = ■ ■ ■ = ^4 , then, we have three extra 
independent relations 

X^ = \^ = X^ = X^ . (19) 
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Combined with equations dHj) and (IT^ . we get A^ = • • • = A^ = 
which is to be interpreted as /^ = ■ ■ ■ = /^ = , see equation (ITB]) . Hence, 
as described earher, it follows for i = l,---,7 that if AJ = initially at 
t = to then Aj = and A* = f* for all t where v'^ are constants. Or, 
if e^ ~ t^ — 7- oo in the limit t — )■ oo with /3 > 1 , it then follows for 
i = 1, ■ ■ ■ , 7 that AJ — 7- and A* — t- v^ in this limit. Obtaining the values 
of the asymptotic constants f * , however, requires knowing the details of 
evolution. It also follows similarly that e^' ~ e s" — >■ oo for i = 8, 9, 10 . It 
is straightforward to show that same results are obtained for the equivalent 
55'5"W configuration also. 

Thus, assuming either that A^ = ■ ■ ■ = Aj = initially at t = to or that 
e^ ~ f'^ — > oo in the limit t — )■ oo with /3 > 1 , we obtain that the A^ = 4 
mutually BPS intersecting brane configurations with identical numbers of 
branes and antibranes, i. e. with rii = ■ ■ ■ = n^ and ni = ■ ■ ■ = n^ , will 
asymptotically lead to an effective (3 + 1) dimensional expanding universe 
with the remaining seven spatial directions reaching constant sizes. This 
result follows as a consequence of U duality symmetries alone, which imply 
relations of the type given in equation ( TT6l) among the components of the 
energy momentum tensor Tab ■ This result is otherwise independent of the 
details of the equations of state, and also of the ansatzes for Tab we make 
in the following in order to proceed further. 

C. Ansatz for Tab 

The dynamics underlying the general result given above may be under- 
stood in more detail, and the asymptotic constants v^ can be obtained, if an 
explicit solution for the evolution is available. In the following, we will make 
a few assumptions which enable us to obtain such details. 

Consider now the case of 2 branes or 5 branes only. From the U duality 
relations given by equations fITT]) . flT^ . flT51) . and f lTB]) . it follows easily that 
P|l = —p + 2pj_ where p\\ is the pressure along the brane directions and 
p± is the pressure along the transverse directions. For the case of waves, 
one obtains p\\ = p ■ We write these U duality relations for the A^ = 1 
configurations in the form 

P\\= z{p-p±)+p± (20) 

where z = — 1 for 2 and 5 branes and = +1 for waves. (A similar relation 
may be obtained in the black hole case also, see Appendix A.) In general, p , 
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p\l , and p± are functions of the numbers n and n of branes and antibranes, 
satisfying the U duahty relations (I2UI) . If n = n then p\\ and p± may be 
thought of as functions of p satisfying equation fl20|) [24J. 

Consider now mutually BPS A^ intersecting brane configuration. In the 
black hole case, it turns out that the energy momentum tensors T^^^js of 
the J*'* type of branes are mutually non interacting and seperately conserved 
Hg - [SI]. That is, 

n = En(/) , EVAn(/) = o . (21) 

I A 

We assume that this is the case in the context of early universe also where 
T^B = diag{-p, pi) , T'^Bii) = diag{-pj, pu) , p/ > , and (p/, pu) satisfy 
the U duality relations in (120|) for all / . Equations (!2T]) now give 

P = J2pi ^ Pi = Y.Pii (22) 

{pi)t + piAt + Y.PiiK = ■ (23) 

i 

We have verified explicitly for a variety of mutually BPS A^ intersecting 
brane configurations that equations fl20l) and fl22|) are sufficient to satisfy all 
the relations of the type J2i q/* = implied by U duality symmetries. See 
[Mj for more details. 

To solve the evolution equations ([7]), ([H]), (122]) . and ( 1251) . we need the 
functions pj , p\\j , and p^/ . To proceed further, we assume that nj = nj 
for all / . This is necessary if, as is the case here, the common transverse 
directions are compact and hence the nett charges must vanish. Then p^j 
and p±i may be thought of as functions of pi satisfying equation (120|) . 

It is natural to expect that p±i{pi) is the same function for waves, 2 
branes, and 5 branes since they can all be transformed into each other by U 
duality operations which do not involve the transverse directions. We assume 
that this is the case. We further assume that this function p±{p) is given by 

p^ = il-u)p (24) 

where m is a constant. Such a parametrisation of the equation of state, 
instead of the usual one p = wp , leads to elegant expressions as will become 
clear in the following, see [321 [33] also. The results of [20] correspond to the 
case where u = 1 . Here, we assume only that < tt < 2 . The constant u is 
arbitrary otherwise. 
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It now follows that pu in equation ([22]) are of the form pn = (1 — wf ) pi 
and that the constants wf can be obtained in terms of u using equations 
f l20|) and fl24|) . Thus, for 2 branes, 5 branes, and waves, we have u± = u , 
u\\ = (1 — z) u , and hence 



2 

5 

W 



Ui 
Ui 



2, 2, 1, 1, 1, 1, 1, 1, 1, 
2, 2, 2, 2, 2, 1, 1, 1, 1, 
0, 1, 1, 1, 1, 1, 1, 1, 1, 



where the / superscripts have been omitted since A^ = 1 
the 22'55' configuration are given by 



1) u 
1) u 
1) u (25) 

Similarly, ul for 



2 

2' 
5 
5' 






2, 
1, 
2, 
1, 



2, 1, 1, 

1, 2, 2, 

1, 2, 1, 

2, 1, 2, 



1, 1, 1, 



1, 
2, 
2, 



1, 
2, 
2, 



1, 1, \)u 

1, 1, l)n 

1, 1, \)u 

1, 1, 1)m . 



(26) 



This completes our ansatz for the energy momentum tensor T^^ for the 
intersecting brane configurations in the early universe. 



III. General Analysis: Evolution equations 

The evolution of the universe can now be analysed. In this section, we first 
present the analysis in a general form which is applicable to a D dimensional 
homogeneous, anisotropic universe. We specialise to the intersecting brane 
configurations in the next section. 

The D dimensional line element ds is given by equation (jH]), now with 
i = l,2,---,D — 1. The total energy momentum tensor T^^ of the dominant 
consituents of the universe is given by equation ([6]). The equations of motion 
for the evolution of the universe is given, in units where SttG^) = 1 , by 
equations ([7]) - ([9]) with 9 in equation (|8]) now replaced by D — 2 . Defining 



Gij 



\-b., 



ij 1 



G'^ 



D 



2 ^'' ' 



(27) 
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the equations ([7]) and ([H]), with 9 replaced by D — 2 , may be written com- 
pactly as 



Y.G.jXlXl = 2p (28) 

i 

Xl, + AAl = Y.G'Hp-Pj) (29) 



^,J 



where i, j,- ■ ■ run from 1 to D — 1 . 

Let the universe be dominated by A^ types of mutually non interacting 
and seperately conserved matter labelled by / = 1, ■ ■ ■ , A^ . Then the corre- 
sponding energy momentum tensors Tab{i) and their components pi and 
Pii satisfy equations f l2T]) - f l23|) . 

Further, let the equations of state be given by pu = (1 — itf) p/ where 
ul are constants. Equations fl25]) . fl25]) . and fl2^ may now be simplified and 
cast in various useful forms as follow. 

Using pii = {1 — ul) Pi , equation ( 123|) can be integrated to give 

P/ = e''-^^ , l' = Y.<>^ + ll (30) 

i 

where /q are integration constants. Further using equations ( 12^ and (pl]l . 
equations (|28|) and ( 129|) become 

Y.G.,K^t = 2Y.e"-'^ (31) 

Kt + ^tK = Y.^''e''-'^ (32) 

where ti*"' = X]j G^^ uj . Let the initial conditions at an initial time t^ be 
given, with no loss of generality, by 

{\\ XI, I', ll pi)^^^^ = (O, k\ ll K\ p,o) (33) 

where 

p,o = e^° , K' = Y^ uik' , Y. G^Jk^k^ = 2 ^ e'o . (34) 

i i,j J 

Equations f l3T|) and f l32|) may now be solved for the D—1 variables A* with the 
above initial conditions. Or, instead, these equations may be manipulated 
so that one needs to solve for A^ variables l^ only, see equations fl35|) . fl39|) . 
and (jH|) below. We now perform these manipulations. 
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First define a variable T{t) as follows: 

dr = e-^ dt , r(to) = . (35) 

Then, for A*(t) or equivalently A*(T(t)) , we have 

XI = e~^ K , Kt + ^tK = e~'^ Kr (36) 

where the subscripts r denote r-derivatives. Note that the initial values at 
r(to) = remain unchanged since A = , and hence A^ = A^. at t = to . 
Equations fl3T|) and fl32|) now become 

Y.G^,KK = 2Ee'' (37) 

Kr = E^''e^' ■ (38) 

J 

Also, from l^ = J2i ujX'' + l^ , it follows that 

C = T.S''e^' (39) 

J 

where 

G'' = E^l^' = EG''<^' ■ (40) 

We assume that Qu exists such that J2.jGij G"^^ = Sj^ , i.e. that the 
matrix Q formed by Q^"^ is invertible. g Then, from equation fl39l) . we have 

T.Sijt = e'' ■ (41) 

J 

Substituting this expression for e' into equation (!38|) . then integrating it 
twice and incorporating the initial conditions given in equation (1331) . we get 

A^ = E «/ (^' - ^o) + LW , u} = Y,gjj G'^ uj (42) 



^This is not always the case. For example, uf — u^ for all i in the isotropic case. 
Then Q^"^ oc u^w^ and Q is not invertible. This is not a problem, it just means that the 
set of variables l^ can be reduced to a smaller independent set; one then proceeds with 
the smaller set. 
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where U are integration constants. It follows from A^O) that U are related 
to initial values /c* and K^ by k^ = J2i u\ K^ + U . Using this expression 
for fc* in the relation K^ = J2i u^k^ , or substituting the expression for A* 
given in equation P2l) into the equation fl3Ql) for I' , leads to the following A^ 
constraints on L* : 

Y,ulV = , I=1,2,---,N . (43) 

i 

Now, using equations (142|) and (HSj) . equation (!37|) may be written in terms 
of F as follows: 

Y: Qij lili = 2{E + Y^e'') , 2E = -Y: G,,UU . (44) 

i,j I ij 

One may now solve equations (139|) and (jH]) for A^ variables 1\t) . Then A*(t) 
are obtained from equation fH21) and t(r) from equation (15^ . Inverting t(r) 
then gives T(t) , and thereby A*(t) . 

A. A^ = 1 case 

Consider the A^ = 1 case. Note that we are still considering the general 
D dimensional universe, not the eleven dimensional one. We assume here 
that Q^^ = Q > . Now, as shown in Appendix B, it follows in general that 
if J2i UiU = and J^ij G'^^UiUj > then E >0 and E vanishes if and only 
if L* all vanish. Since Y.i^lL'^ = , see equation fH51) . and we assume that 
Q^^ = J2i,j G'^uju] > , we have E > . We further assume that E > , 
equivalently that L*s do not all vanish. 

Omitting the / labels, equations (139!) and (jH]) for /(r) become 

l,, = ge' , {lrf = 2g{E + e') . (45) 

The initial values are k = /(O) and K = /^(O) obeying K"^ = 2 g {E + e'") . 
We take K > with no loss of generality. Then the solution for /(r) is 
given by 

Sinh^ a{Too — r) 
where 

2a2 = gE , Sinh^ ar^ = E e"'" , K = 2a Coth ar^ . (47) 
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The sign of a is immaterial but, just to be definite, we take it to be positive. 
The sign of t^q is same as that of K , hence Tqo > . Also, A'(r) and ^(t) 
may now be obtained but are not needed here for our purposes. 

Note that e' — )■ 4£' e^"*^'^~'^°°) and vanishes in the limit r — )■ — oo , 
whereas e' — > ^ (tqo — r)^^ and diverges in the limit r — )■ Too . The value of 
Too depends on the initial values ^o and K, or equivalently E , as given in 
equations fll7|) . It is finite and can be evaluated exactly. However, if e'° ^ E 
then Too niay be approximated in a way that will be useful later on. 

From the exact solution given above, we have Sinh"^ ar^o = E e"'" and 
K = 2a Coth oToo . In the limit e'-° <t: E , we then have e^"^°° ~ 4 £" e~'° 
and K c:^2a . It, therefore, follows that 

Too ^ —{lnE-k + In 4) . (48) 

K 

In the limit e''° <t^ E , the evolution of /(r) may also be thought of as 
follows. Consider E to be fixed and e'° to be very small so that e^° <^ E . 
It then follows from equations ( 145|) that, at initial times, Irr is very small 
and that Ir — \/2QE = 2a is independent of e' . Hence, 1{t) evolves as if 
there is no 'force', i.e. 1{t) ~ /q + Kt where K = It{0) > is the initial 
'velocity'. Once e^ becomes of 0{E) then it affects Ir ■ But, from then on, 
e' evolves quickly and diverges soon after. 

This suggests that one may well approximate Too by the time Tq required 
for / , which starts from ^o with a velocity K and evolves freely with no 
force, to reach In E - namely, to reach a value where e' = 6^°^^"^" = E . In 
other words, if e'" -C E then 

T^c:^Ta= -^7 {In E - Iq) . (49) 

A 

A comparison with equation fj48|) shows that the exact t^o which follows from 
solving the evolution equations is indeed well approximated by r^ in equation 
( 149|) in the limit e^° '^ E . Note that r^ is calculated using only the initial 
values, requiring no knowledge of the exact solution. 

B. A^ > 1 case 

When A > 1 , the equations of motion can be solved if Q^^ are of certain 
form [32] - |38]. For example, if Q^'^ oc 5^^ then the solutions are similar 
to those in the A^ = 1 case described above. For general forms of Q^^ , we 
are unable to obtain explicit solutions. Nevertheless, the general evolution 
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can still be analysed if one assumes suitable asymptotic forms for the scale 
factors e^* . 

It follows from equations f l27|) and f l28|) that At cannot vanish. With 
no loss of generality, let Aj > initially at t = to • Then e^ decreases 
monotonically for t < to , equivalently r < , and increases monotonically 
for t > to , equivalently r > . Further features of the evolution depend on 
the structure of Q^"^ and uj . In the cases of interest here, it turns out that 
e^ and also all e' vanish in the limit r — )■ — oo , and diverge in the limit 
'T ^ Too where r^o is finite. We assume such a behaviour and analyse the 
asymptotic solutions. 

1. Asymptotic evolution: e^ —^ 

We assume that (e^, e' ) — )■ in the limit r — )■ — oo . Then, equations 
f pSjl and fl59]) can be solved since their right hand sides depend only on e' s 
which now vanish. Hence, in the limit r — )■ — oo , we write 



-') 



which are valid upto multiplcative constants and where (c^, c*, b^ , U) are 
constants. Also, the equalities in the asymptotic expressions here and in the 
following are valid only upto the leading order. Equation fH2]) now implies 
that c* = X]/ W/C^ + L^ ■ Also, e^ = e^"^ where c = X]j c* . Then it follows 
from equation (135!) that t ~ e*^"^ . Hence, 



lI 



C ~i c 



c c 



b' = ^ , T.b' = l ■ (51) 



Furthermore, equation (1571) implies that (I]j&*)^ — Z]j(^*)^ = . Thus the 
evolution is of Kasner type in the limit r — t- — oo . The constants c^s in 
equations fl50|) must be such that the resulting J2i ^* = J2i{b^Y = 1 ? but are 
otherwise arbitrary. In an actual evolution, however, c^s can be determined 
in terms of the initial values /q and K' with no arbitrariness, but this 
requires complete solution for 1\t) . 

2. Asymptotic evolution: e^ — ?■ oo 

We assume that e^ — t- oo in the limit r — )■ Too where Tqo is finite. 
Whether this limit is reached at a finite or infinite physical time t depends 
on the values of uf , see below. A(r) may be obtained in terms of 1^{t) using 
equation fH2]) . Hence, in the limit e^ — > oo , some or all of the e' s diverge. 
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Consider the following ansatz in the limit r — )■ Tqo : 

e'' = e^' {r^ - r)-^'' , e^" = e^" {r^ - r)~'^' , (52) 

where c^ and 7^ are constants, and some or all of the 7''^s must be > so 
that some or all of the e' s diverge. Equation ( l42l) now implies that 

f = E«7 7' , c^ = E«/ (c' - ^0) + LW^ . (53) 

Also, e^ = e'^ {toq — t)^^'^ where c = J^iC' and 7 = Z)i 7* • For the ansatz in 
equations fl52|) to be consistent, it is necessary that 7 > so that e'^ — )■ 00 
in the limit r — ?■ Too • Now ^(t) follows from equation ([3SD and is given by 

t - ts = ^-^ e^ {t^ - rp'^''^ , 7 = E«1/ (54) 

where t^ is a finite constant. If 27 < 1 then t ^ tg which means that 
e^ — )■ cxD at a finite physical time tg . If 27 > 1 then t — )■ 00 in the limit 
e^ —7- 00. Which case is realised, i.e. whether 27 < 1 or > 1 , depends on 
the values of ul . 

Using equation ( 15^ . the asymptotic behaviour of e' and e^' can be 
obtained in terms of t. For example, let 27 > 1 and 

/ = e''^^' t^' , e^' = e"" t^' , e^ = e' t^ (55) 

in the limit t — > 00 . It then follows that 

,3'=^^. ^'=^^, /3=^^. (56) 

^ 27 - 1 ' ^ 27 - 1 ' ^ 27-1 ^ ' 

Note that, in this case, we have e^ ~ t^ in the limit t — )■ 00 with /3 > 1 . 
See the discussion below equation (iT6l) for the relevance of this feature. 

To obtain the values of 7^ , and thereby 7* , in equation fl52|) . consider 
equation fj4T]) from which it follows that 

2Y.QiJl'= e'^' {r^ - rf^'-''^ . (57) 

The left hand side in the above equation is a constant but the right hand 
side depends on r . This is consistent if 7^ = 1 in which case the right hand 

21 



side becomes a positive constant, or if 7"^ < 1 in which case the right hand 
side vanishes in the hmit r — )■ Tqo • Thus, there are two possibihties: 

(i) y = i ^ 2Y.gul' = e^' >Q (58) 



J 
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t't^I ^ Y.'3iJl' = ^ , 7'<1 • (59) 



J 



For a given Qu , the possible consistent solutions for (7^, e'^ ) are to be 
obtained as follows. Assume that some J's are of type (i) and the remaining 
ones are of type (ii) . Then solve equations fl58l) and fl59|) for e^ in type (i) 
and for 7^ in type (ii). Such a solution is consistent if the resulting (e*^ , 7^) 
satisfy e'^ > for /s in type (i) and 7^ < 1 for Is in type (ii). Also, some 
or all of the 7^s must be > as required in equation ([S2D- (It is further 
necessary that the resulting 7 > so that e"^ — )■ 00 , but calculating 7 
requires ul .) 

Consider an example, which will be useful later, where Q^^ and Qu are 
given by 

G^' = a{b~6^') , Gij = -{-^-5u) (60) 

a i\b — 1 

with a > and Nb > 1 . It is then easy to show that the only possibility 

is the one given in (i). Also J2jQiJ = a(Nb-i) ^ ^ ' ^^^ ^'^^^ 7^ = 1 for 

all J is a consistent solution as required by equation ( l58l) . In the A^ = 1 

case, we get Q^^ = Q = a{b— 1)>0, and e' in the limit r — > Too obtained 

as described above agrees with that obtained from the explicit solution, see 

below equation (H71) . 

Thus e'^ and 7^, and thereby 7* = J2i u) Y and 7 = J2i,i u} 7^ , are all 

determined ultimately by uj . The constants c* are given by equation (l53l) 

and they depend on u^ , on the initial values l^ and L* , and also on Too ■ 

But determining t^o , and hence determining c* when L* do not all vanish, 

requires complete solution for 1\t) . 

C. Deviations from e' — )■ 00 Asymptotics 

We consider the deviations of 1\t) from its asymptotic behaviour given 
in equation fl52|) . which will turn out to be of interest. Let the deviations 
s^{t) ioT I = 1,2,- ■ ■ ,N be defined, in the limit r — )■ Too , by 

e'' = e^' (Too - r)-'''' e^^^ (61) 
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where c^ and 7^ are determined as described earlier. For the purpose of 
illustration, and also for later use, we now assume that all the /s are of type 
(i), namely that 7^ = 1 and e'^ = 2 Yl,j Qu > for all / . It then follows 
straightforwardly from equation fl39|) that 



r^ 



rf 4. = 2 E G'"" Qkl (e^" - 1) • (62) 



K,L 



Consider the example of Q^"^ given in equation (160|) . Then J2jGij 
^ and, for any a^ , one has 



a(7V6-l) 



E^'^^^.-^ = -T^(-'-^E-^) ■ (63) 

K,L IMU I ^ 

In equation fl62l) . a^ = 2 (c' ~ 1) • I^ ^ow follows easily that, up to the 
leading order in {s^} , the difference s^ — s"^ obeys the equation 

(roo - rf (/ - s')^^ + j^^ (/ - .^) = . (64) 

The solutions to these equations are of the form 

(^^-^^) ~ (^oo-r)^(^^v^) , A = l-^^. (65) 

Note that s^ — s'^ = l^ — F since 7^ and c^ are same for all I , see equation 
fl6T|) . Hence, equations fl64l) and fl65|) can be used to understand in more 
detail the behaviour of l^s as they all diverge in the limit t ^ Too as given 
in equation fl52|) . We will discuss these features in sections VI and VII. 



IV. Intersecting Branes 

We now analyse the evolution of the universe dominated by mutually BPS 
A^ intersecting brane configurations of M theory. The number of spacetime 
dimensions D = 11 . The equations of state are assumed to be given by 
Pii = (1 — uj) Pi where, as a consequence of U duality symmetries, u^ are 
parametrised in terms of one constant u . The indices i,j, ■ ■ ■ run from 1 to 
10 and the indices I, J,- ■ ■ from 1 to A^ . For 2 branes, 5 branes, and waves, 
A^ = 1 and the corresponding ul are given in equations fl25|) . For 22'55' 
configuration, A^ = 4 and the corresponding u^ are given in equations 
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A. Evolution Equations 

The evolution of A* describing the scale factors is given by the equations 
described earlier which, for ease of reference, we summarise below: 







Kr = 


J 


(66) 






lir = 


J 


(67) 






y = 


^ U^j {I' - li) + V T 
J 


(68) 


ire 




, G'' = 




(69) 



with G"^^ and Qu as defined earlier, and U are arbitrary constants satisfying 
the constraints X^i '^{L^ = for all / . Also, /^ obey the constraint 

YGijlili = nE + Ye'') (70) 

i,j J 

where 2E = — J^ij GijUU . Equations (!67j) and flTOj) are to be solved for 
1\t) with initial conditions /^(O) = Iq = In p/o and /:^(0) = K^ where pio 
are initial densities and 

YGijK' K' = 2{E + Ye''') ■ (71) 

i,j J 

Then A*(t) follow from equation (!68|) and the physical time t(r) from 
dt = e^ dr . Inverting t^r) then gives r(t) , and thereby A*(t) . 

We can now calculate Q^'^ for the mutually BPS intersecting brane con- 
figurations. As explained in footnote 2, in the BPS configurations two stacks 
of 2 branes intersect at a point; two stacks of 5 branes intersect along three 
common spatial directions; a stack of 2 branes intersects a stack of 5 branes 
along one common spatial direction; and, waves, if present, will be along a 
common intersection direction. With these rules given, it is now straightfor- 
ward to calculate Q^'^ using equations ( l25l) and ( l69l) . It turns out because 
of the BPS intersection rules that the resulting Q^'^ are given by 

g'^ = 2u^ (1 - 5'^) . (72) 
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The corresponding Qjj exists for A^ > 1 , and is given by 

1 



Gij 



1 
2^ 



•A^ 



6 



iJ 



(73) 



Note that, for A^ > 1 , the above Q^'^ is a special case of the example 
considered earlier in equation fl60|) . now with a = 2v? and 6=1, 

It is also straightforward to calculate li*^ and u\ for the 22'55' config- 
uration using the definitions in equation flB^ and the u\ in equation 
They are given by 



2 


w^i 


oc 


(-2, 


-2, 


1, 


1, 


1, 


1, 


1, 


1, 


1, 


1) 


2' 


^'' 


oc 


( 1, 


1,- 


-2, 


-2, 


1, 


1, 


1, 


1, 


1, 


1) 


5 


u'' 


oc 


(-1, 


2, 


-1, 


2, 


-1, 


-1, 


-1, 


2, 


2, 


2) 


5' 


u^' 


ex 


( 2, 


-1, 


2, 


-1, 


-1, 


-1, 


-1, 


2, 


2, 


2) 



(74) 



where the proportionality constant is ^ , and by 



2 
2' 
5 
5' 



ni oc ( 2, 2,-1,-1,-1,-1,-1, 
n^2 oc (-1,-1, 2, 2,-1,-1,-1, 



n, oc 



-2, 1,-2, 1, 1, 1, 



1, 1, 1) 
1, 1, 1) 
0, 0, 0) 



< oc (-2, 1,-2, 1, 1, 1, 1, 0, 0, 0) 



(75) 



where the proportionality constant is ^ . 

We are unable to solve equations f E7|) . (170]) . and (172]) for A^ > 1 o How- 
ever, applying the general analysis described in section III and making fur- 
ther use of the explicit forms of ttf and Q^'^ given in equations (126]) and (172]) . 
one can understand the qualitative features of the evolution of the 22'55' 
configuration. 

We first make several remarks which will lead to an immediate under- 
standing of the evolution of this configuration. 



^ In the case of black holes, the equations of motion for the corresponding harmonic 
functions H^ = 1 + — = e^' can also be written in a form similar to that of equation 
(|67| . The main steps are indicated in Appendix A. The analogous Q^'' in the black hole 
case turns out to be oc 5 , and the equations can then be solved. 

Also, note that if L* = for all i then A* in equation (|55)) here may be written 
as in equation (|108p in Appendix A. The role of hj there is played by the functions 
2uhj — 2u ^ J Gij {l'^ — Iq) here. Such a similarity is present for other intersecting brane 
configurations also. 



25 



(1) Let Ui = J2iul ■ It can then be checked that Y.i,jG'^^UiUj > . 
Also, X^i i^iL^ = since J2i ^f ^* = for all / . Hence, as shown in Appendix 
B, it follows that E given in equation flTOj) is > and that it vanishes if 
and only if L* all vanish. 

(2) The constraints X^i "wf-^* = imply that 

^8 + L^ + L^o = -3 (L^ + L^) . (76) 

Thus, for example, we may take {L^ , L'^ , L^ , L"^ , L^ , L^) to be independent. 
The remaining L*s are then determined by the above equations. Also, we 
have 

L = ^L' = -(L^ + L2) . (77) 

i 

Using equations (175]) . (177|) . and the Schwarz inequality (11201) in Appendix 
B, we write E as 

2E = E(^T-(E^T 

i i 

= ?>{L^f + {L^ + 2L''f + Y,{Uf + al (78) 

where the first line is the definition of i? , cr2 = if and only if L^ = L"^ , and 
(Ts = if and only if L^ = L^ = L^° . See the Schwarz inequality given in 
equation (11201) . It is easy to show that the above expressions for E imply that 
(L*)^ for all i are bounded above by E as follows: E > Ci{U'Y — where 
Cj are constants of C(l) . In particular, note the inequality 2E > 3(L)^ 
which is required in Appendix C. 

(3) It follows from equations (!68l) . (175|) . and (!77j) that 

A- = E-^^ = ^(2C + 2/? + ^' + a + ^ ■ (79) 



Using the explicit form of Qjj given in equation (1731) with A^ = 4 , equation 
rol) becomes 

(E ^ir - 3 Ei^ir = 12t.^ iE + j:e'') > O (so) 
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where the inequahty follows since E > and e' > . We show in Appendix 
C that this inequality implies that none of (A^, l^^) may vanish, and that 
they must all have same sign. Hence, for all r throughout the evolution, 
(A-r, l^) must all be non vanishing, and be all positive or all negative. 

B. Asymptotic evolution 

With no loss of generality, let A^ > initially aX t = to . Then it follows 
from the above result that (A,-, l^) must all be positive and non vanishing 
for all T . Hence, (A, l^) are all monotonically increasing functions for all r 
throughout the evolution. 

Equation fl67|) may be written, using equation f l72|) . as 



TT 



2u^ Y, e' . (81) 



In the past, r and all l^ decrease continuously. Hence, the right hand side 
in equation ( IHTj) becomes more and more negligible. The asymptotic solution 
in the limit r — t- — oo is then given by l^ = cW + S where c^ > . Thus 
e' — )■ in this limit. 

Similarly, in the future, r and all f increase continuously. However, 
the right hand side in equation (IHTj) increases exponentially now. It is then 
obvious that all e' — )■ cxd within a finite interval of t , i.e. at a finite value 
Too of r . In this context, see equations (H5!) and (146!) . and the general 
analysis given in section III B.2. 

We now analyse the corresponding asymptotic solutions. 

1. Asymptotic evolution: e"^ — > 

It follows from the above discussion that e^ — t- in the limit r — > — cxd . 
Also, in this limit, we have 



ii .--f. 



,b' „X^ ^C^T J) 



e'^ = e" " = r , e^ = e"^ = f (82) 

upto multiplcative constants where (c'^, c\ h^ , U) are constants. The 
evolution is then of Kasner type and is similar to that described in section 
III B.l. The constants c'^s are determined by the initial values /q and K^ , 
but obtaining the exact dependence in the general case requires complete 
solution for 1\t) . However, if the initial values /q are large and negative 
then we have e' ^ 1 for all r < and, hence, c^ = K^ to a good 
approximation. 
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2. Asymptotic evolution: e^ — )■ oo 

It follows from the above discussion that e^ — > cxd in the limit r — )■ Too 
where t^q is finite. Also, e' — > oo in this limit and t^q depends on the 
initial values /q and K^ . 

Although solutions for 1\t) are not known, their asymptotic forms in 
the limit r — )■ Too , and hence those of A*(t) , may be obtained following the 
analysis given in section III B.2. Q^'^ in equation (I72|) is a special case of 
the example fl60l) where, now, N = A , a = 1v? , and h =\ . Hence, it can be 
shown to correspond to the possibility (i) given in equation (1551) . Therefore, 
we have 7-^" = 1 and e''' = 1Y,jQij = ^ ■ 

It then follows from equation ( 152|) that e^ and e'^' are given in the limit 
T ^ Too by 

e" = ^.7-^. (83) 



e^' 



3^2 (Too - tY 

1 \E/"} 



(84) 



where, since p/o = e'o , we have 

v' = -Y.^'jli + UToo , e^' = e'^'^- n(PJo)-"'^ • (85) 

J J 

Also, since 7 = J2i,iu\ = ^ , we have from equation ( 15^ that the physical 
time t is given in this limit by 

2—11 

t-ts = A (Too- r)- — (86) 

where tg and A are finite constants. Clearly, t — ;■ 00 in the limit t —^ Too 
since it is assumed that < m < 2 . In this limit, the scale factors e^' may 
be written in terms of t as 



e 



A' 



e"" (Bt)^' (87) 



where _B is a constant and /?* = ^^ X) j "^j • Using equation (175|) for «} , the 
exponents /3* are given by 

/3' oc (0, 0, 0, 0, 0, 0, 0, 1, 1, 1) (88) 
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where the proportionahty constant is g,^^ s . Note that /3 = I]i /3* = 23^ > 
1 . Hence, we have e^ ~ t^ in the hmit t — t- oo with /3 > 1 . See the 
discussion below equation ( IT6|) for the relevance of this feature. 

• i 2 

Thus, asymptotically in the limit t — )• oo , we obtain that e — )■ t3(2-u) 
for the common transverse directions i = 8, 9, 10 . Hence, these directions 
continue to expand, their expansion being precisely that of a (3 + 1) - di- 
mensional homogeneous, isotropic universe containing a perfect fluid whose 
equation of state is p = {1 — u) p . Also, e^' — )■ e^' for the brane directions 
i = 1,- ■ ■ ,7 . Hence, these directions cease to expand or contract. Their sizes 
are thus stabilised and are given by e'"' . Note that this result is in accord 
with the general result described in section II B since, in the limit r — )■ Tqc, , 
the brane densities pi oc e' all become equal and hence the four types of 
branes all become identical; and, t — )■ oo and e^ ~ t^ — )■ oo with (3 > 1 . 

C. Mechanism of stabilisation 

Using the asymptotic solutions, we can now give a physical interpreta- 
tion of the dynamical mechanism underlying the stabilisation of the brane 
directions seen above for the 22'55' configuration. 

We first study the stabilisation process. Consider equation (166|) for A^^ , 
for example. Using the values of u^^ given in equation (174|) . we have 

A^, oc (-2e'' + e'" - e'" + 2e'*) . (89) 

In the 22'55' configuration, x^ direction is wrapped by 2 branes and 5 branes 
and is transverse to 2' branes and 5' branes. Thus, from the above equation 
for A^ and from similar equations for A^, ■ ■ ■ , A^ , we see that 2 brane and 
5 brane directions 'contract with a force' proportional to 2p(2) and p(5) 
respectively, whereas the directions transverse to them 'expand with a force' 
proportional to p(2) and 2p(5) respectively, where p(*) oc e' * are the time 
dependent densities of the corresponding branes. 

When pj oc e' all become equal, the forces of expansion cancel the forces 
of contraction resulting in vanishing nett force for the x^ direction. Then, 
using equation (15B]) . one has 

Xl^ = e'^ {XI, + AA]) = ■ (90) 

Now, as described earlier in the context of equations (IT7|) and (I18p . the 
transient 'velocity' A^ is damped and A^ reaches a constant value in the 
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expanding universe here since we have e^ ~ t^ in the hmit t — t- oo with 
P > 1 . The result is the stabihsation of the x^ direction. 
The stabihsed size e^ of x^ direction is given by 

e^'=e^'^- ((^f , (91) 

VPsoPio/ 

see equation (!85l) . Note that e"" can be interpreted as arising from the 
imbalance among the initial brane densities pio , and from the parts L^ of 
A^(0) which indicate the transients. The above analysis can be similarly 
applied to the stabilisation of other brane directions (x^,---,x^) in the 
22'55' configuration. 

Thus, three conditions need to be satisfied for stabilisation: (1) the time 
dependent brane densities pi oc e' all become equal; (2) the forces of 
expansion and contraction for each of the brane directions be just right so 
that the nett force vanishes; (3) the universe be expanding as e^ ~ t^ in 
the limit t — > oo with /3 > 1 so that the transient velocities are damped 
and the corresponding scale factors reach constant values. 

For any mutually BPS A^ > 1 intersecting brane configurations with 
the equations of state as assumed here, it is straightforward to show using 
the earlier analysis that the evolution equations ensure that e' all become 
equal asymptotically even if they were unequal initally, and that e^ ~ t^ in 
the limit t — )■ oo with /3 > 1 . Thus conditions (1) and (3) are satisfied. 
Condition (2) requires the brane configuration to be such that each of the 
brane directions is wrapped by, and is transverse to, just the right number 
and kind of branes. This condition is satisfied for the A^ = 4 configurations 
22'55' and 55'5"W , both of which result in the stabilisation of seven brane 
directions and the expansion of the remaining three spatial directions. To 
our knowledge, the only other configurations which satisfy the condition (2) 
are the A^ = 3 configurations 22'2" and 25W , both of which result in 
the stabilisation of six brane directions and the expansion of the remaining 
four spatial directions [21]. However it is the A^ = 4 configurations that are 
entropically favourable, see equation ([T]). 

Note that, as described in section II B and upto certain technical assump- 
tions regarding the equality of brane densities and the asymptotic behaviour 
of e^ , the stabilisation here follows essentially as a consequence of U duality 
symmetries. In particular, it is independent of the ansatz for energy momen- 
tum tensors, or of the assumptions about equations of state, as long as the 
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components of the energy momentum tensors obey the U duahty constraints 
of the type given in equation (IT^ . Obtaining the details of the stabihsation, 
however, requires further assumptions e.g. of the type made here. 

Note also that the present mechanism of stabilisation of seven brane di- 
rections, and the consequent emergence of three large spatial directions, is 
very different from the ones proposed in string theory or in brane gas models 



V. Stabilised sizes of brane directions 

We thus see for the 22'55' configuration that, asymptotically in the limit 
e^ — 7- oo , the initial (10 + 1) - dimensional universe effectively becomes 
(3 + 1) - dimensional. Also, if v^ = min{v^, ■ ■ ■ ,f^} then a dimensional 
reduction of the (10 + 1) - dimensional M theory along the corresponding x* 
direction gives type llA string theory with its dilaton now stabilised. Using 
the standard relations, one can obtain the string coupling constant gs , the 
string scale Ms , and the four dimensional Planck scale M4 in terms of the 
M theory scale Mn and the stabilsed values e'"' . Defining v'^ = Y.J=i v'^ and 
assuming, with no loss of generality, that the coordinate sizes of all spatial 
directions are of 0{M{i) , we obtain 

^2 = e3-^ , M^ = e^'-^' M^ = e^' M^, (92) 

where the equalities are valid upto numerical factors of (9(1) only and 

VP3oP4oy ' h 

as follows from equations (^^, fl55]) . and pjo = e'o . Also, note that gs = 

( Ms \3 
V A/11 ^ ' 

Since we have an asymptotically 3 + 1 dimensional universe evolving from 
a 10 + 1 dimensional one, it is of interest to study the resulting ratios ^j^ 
and ^ , and study their dependence on the initial values {ll,K^,D) . In 
particular, one may like to know the generic values of these ratios and to 
know whether arbitrarily small values are possible. Setting M4 = 10^^ GeV , 
one then knows the generic scales of Mn and Ms and, for example, whether 
Mil = 10~^^ M4 = 10 TeV is possible. 

In view of the relations between {Mii,Ms,M4) given in equation (|92|) . 
this requires studying the stabilised values e"" and e^"'""" , their dependence 
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on {Iq,K^ ,U) , and knowing whether they can be arbitrarily large. Note 
that if L* = for all i then f * are all determined in terms of Iq only, 
see equation fl85|) . It is then obvious from equations f l85|) and fl93|) that any 
values for e"" and e^"'""" , no matter how large, may be obtained by fine 
tuning pio correspondingly. [^ 

This statement remains true even when Us do not all vanish. In this case, 
however, one may question the necessity of fine tuning since, for example, 
the relation e"" oc 6'^'"'^°° suggests that large values such as 10^° ~ e™ may 
be obtained by tuning Us, or Too , or both to within a couple of orders of 
magnitude only. It turns out, as we explain below, that fine tuning is still 
necessary to obtain such large values. 

Consider first the possibility of tuning L* . Note that equations (1^ and 
rO|) are invariant under the scaling 



{E, e'^ r) ^ {a'E, a'/, -) (94) 

a 

where o" is a positive constant. The initial values scale correspondingly as 

(e'o, K', V) —^ (a'e'o, aK' , aV) . (95) 

It then follows from equation fl68l) that A*, and hence e""^ , remain invariant o 
This scaling property merely reflects the choice of a scale for time. For 
example, using this scaling, one may set X) j e'" = 1 or, when i? > as is 
the case here, set E = 1 . The corresponding a then provides a natural time 
scale for evolution. We set i? = 1 using the above scaling. 

With E = 1 , the value of Too now depends only on /q and K^ . Since 
2E = J2i{L^y ~ (J2iL^y ; it is still plausible to have a range of non zero 
measure where L* are large and E = 1 , and thereby obtain large values 
for e^" and e"""'""" . However, L*s are further constrained by J2i ^f -^* = , 
/ = 1, ■ ■ ■ , 4 , and consequently their magnitudes are all bounded from above. 
For example, with E' = 1 , we obtain {UY < | . See remark (2) in section 
IV A. Thus, large values of e"* cannot be obtained by tuning L* alone. 



^° It follows from equation ((7T|) and the definition of E that the generic ranges of the 
initial values may be taken to be given by |L*| ~ K ~ ^/E ~ \J pin within a couple of 
orders of magnitude. If the initial values lie way beyond such a range then we consider it 
as fine tuning. 

^^ This invariance is equivalent to that of equations (pS)) and (1^^!) under the scaling 
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Consider now the possibility of tuning r^o ■ Obtaining the dependence of 
Too on {Iq,K^) requires exphcit solutions which are not available. Hence, 
we obtain Tqo numerically. We will present the numerical results in the 
next section. Here we point out that an approximate expression for Too 
can be given in the limit when e'o <^ E for all / . The reasoning involved 
is analogous to that used in obtaining Tq in equation fH^ . Using similar 
reasoning and setting E = 1 now, we have that if e'o ^ 1 for all / then 

Too^Ta= min {tj} , Tj = - ^ . (96) 

Note that r^ can be calculated easily and requires no knowledge of explicit 
solutions. Our numerical results show that Tq given above indeed provides 
a good approximation to Too when e'o ^ 1 for all I . 

Note also that K^ must satisfy equation (1711) with E = 1 . It then 
follows from an analysis similar to that given in Appendix C that K^ are 
all positive, cannot be too small, and are of 0{1) generically. Hence, in the 
limit e'o ^ 1 for all J, r^ in equation (1961) are of (9(mm{— /q}) . This 
indicates that large values of Too , and hence of e"' , cannot be obtained by 
tuning K^ alone; a tuning of /q , which translates to fine tuning of p/o = e'o , 
is required. Our numerical analysis also supports this conclusion. 

We thus find that, even when L*s do not all vanish, a fine tuning of 
P/o = e'o is necessary to obtain large values for e"""'""' and e^" . 

VI. Time varying Newton's constant 

The evolution of the eleven dimensional early universe which is dominated 
by the 22'55' configuration described here can also be considered from the 
perspective of four dimensional spacetime. Indeed, in general, let the eleven 
dimensional line element ds be given by 

7 

ds^ = g^, dx^'dx" + ^ e'^" {dx'f (97) 

i=l 

where x'^ = {x'^ , x^ , x^ , x^^) , with x° = t, describes the four dimensional 
spacetime, and the fields (/^j, and A* , z = 1, ■ ■ ■ , 7, depend on x^ only. Also, 
let A*^ = Yh=i ^* • It is then straightforward to show that the gravitational 
part of the eleven dimensional action 5*11 given in equation ([3]) becomes 

54 = Y^ /rf^xyC^e^^{i?(4) + (V(4)A'=)2-X:(V(4)V)'} (98) 
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where Vy is the coordinate volume of the seven dimensional space and the 
subscripts (4) indicate that the corresponding quantities are with respect 
to the four dimensional metric gfj_i, . The action >S4 describes four dimen- 
sional spacetime in which the effective Newton's constant G4 is spacetime 
dependent and is given by 

G.ixn = e-^^(^') ^ . (99) 

In the case of early universe, the fields g^i, and A* depend on t only. Then 
G4 is time dependent and we have, for G4 and its fractional time derivative, 

G,{t) = e-^^W ^ , ^ = - A^ (100) 

1/7 (jr4 

For the four dimensional spacetime arising from the 22'55' configuration, 
the gfj,,y fields are just the scale factors (e'^ ,e'^,e^ ) for (a;*,x^,x^°) direc- 
tions, and all A*(r) are given in equation (1^ in terms of /^(r) , / = 1, ■ ■ ■ , 4 . 
Then, using equation fl75p and the definitions of A*^ , L'^ , and v'^ , we have 

A-^ = - -i(/l + /2 - /3 _ /4) _ ^c(^^ -t)+v'' . (101) 

DM 

In the limit t — )• 00 , we have from the results given earlier that t ^ Too and 
the fields l^ all become equal. Then A*^ — )■ v'^ where e'"" is given in equation 
(!93|) . and the three dimensional scale factors evolve as in the standard FRW 
case, namely e^ = e^ = e^ ~ ^^(2-") as given in equations fl87|) and fl88l) . 

It thus follows that the effective Newton's constant G4 varies with time 
in the early universe and, in the case of 22'55' configuration, approaches a 
constant value = e~'"'' ^ as the four dimensional universe expands to large 
size. The precise time dependence of G4 will follow from explicit solutions 
to equations (l67|l and (!70|) . The consequences of a such a time dependent G4 
are clearly interesting, and are likely to be important too. But their study is 
beyond the scope of the present paper. 

However, we like to point out here a characteristic feature of the time 
dependence of G4 which arises in the case of 22'55' configuration. Consider 
the behaviour of the differences l^ — l^ in the limit r — )■ Too which, in our 
case, vanish to the leading order. These quantities have been analysed in 
section III C and, for the example of the Q^^ given in equation ( l60l) . they 
are given by equations ( 1M|) and ( l65l) to the non trivial leading order. The 
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case of 22'55' configuration corresponds to A^ = 4 , a = 2u^ , and 6=1. 
Noting that s^ — s'' = l^ — l^ and that A < in our case, equation (1651) 
now gives 

{I'-l^) ~ (Too -r) 5(1 ±^71) (102) 

to the leading order. Clearly, A'^(t) given in equation (1 10 II) will also have 
the same form as above to the non trivial leading order. Thus, taking the 
real part and writing in terms of t using equation (1861) . we have 



A<= = ^;^ + — Sin{u Int + cf)) (103) 

to the non trivial leading order in the limit t — > oo where h and (f) are 

constants, a = 2(2-u) ' ^^"^ ^ ~ v I 2(2-u) ■ Correspondingly, the time 
varying Newton's constant is given by 

G4 oc e"^' = e-"^ (1 - -7^ Sin{u lnt + (j))) (104) 

to the leading order in the limit t — )■ 00 . Note that the constants b and 
depend on the details of matching. The constants a and u arise as real and 
imaginary parts of an exponent on time variable, see equation (11021) . They 
do not depend on the initial values (/g, K^ , U) and thus are independent of 
the details of evolution, but depend only on the configuration parameters A^ 
and u . 

The amplitude of time variation of G4 is dictated by a , and it vanishes 
in the limit t — > 00 . Hence, the time variation of 6*4 in equation (I104p is 
unlikely to contradict any late time observations. The time variation of G^ 
has log periodic oscillations also: G4 has an oscillatory behaviour where the 
77,*'* and {n + l)*'^ nodes occur at times t„ and tn+i which are related by 
In tn+i = ^ + In tn , i-e. by t„+i = e~ t„ . The characteristic signatures 
and observational consequences of such log periodic variations of G4 are not 
clear to us. 

Log periodic behaviour occurs in many physical systems with 'discrete self 
similarity' or 'discrete scale symmetry': for example, in quantum mechani- 
cal systems with strongly attractive -^ potentials near zero energy ^9]; in 
Choptuik scaling and brane - black hole merger transitions |10]; and in a 
variety of dynamical systems [Hj. Algebraically, the log periodicity arises 
when an exponent on an independent variable becomes complex for certain 
values of system parameters. The relevant equations and solutions can often 
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be cast in a form given in equations fl64p and fl65p . But we are not aware of 
a physical reason which explains the ubiquity of the log periodicity. 

To our knowledge, this is the first time a log periodic behaviour appears 
in a cosmo logical context. One expects such a behaviour to leave some novel 
imprint in the universe. But it is not clear to us which effects to look for, or 
which observables are sensitive to the log periodic variations of G4 . 

VII. Numerical results 

We are unable to solve explicitly the equations ( E7}i - ( 170]) describing 
the early universe evolution. Hence, we have analysed these equations nu- 
merically. In this section, we briefly describe our procedure and present a 
few illustrative results. We have analysed both the "U = | and u = 1 cases 
which would correspond to four dimensional universe dominated by radiation 
and pressureless dust respectively. The results are qualitatively the same and, 
hence, we take m = | in the following. Note that u in equation fll04p is then 
determined and, for it = | , the n*^ and (n + 1)*^ nodes in the log periodic 

oscillations occur at times t„ and tn+i related by /n(^y^) = 47rw| ~ 9.734 . 

We proceed as follows. We start at an initial time r = and choose a set 
of initial values /q = In pjq . For each set of /g , we further choose numerous 
arbitrary sets of {K^,U) such that K^ > , E = 1 , and equations ( 17T]) 
and f lTO]) are satisfied, o For each set of initial values (/q, K^ , U) , we then 
numerically analyse the evolution for r > and obtain the value of Too ; the 
evolution of l^ , (A^, ■ ■ ■ , A^°), and t ; the stabihsed values {v^ ,■ ■ ■ ,v'^) ; and 
the resulting values for (^f^, ^^, ||^) . For a few sets of initial values, we have 
analysed the evolution for r < also. 

We find that the numerical results we have obtained confirm the asymp- 
totic features described in this paper: 

(1) e^* and l^ all vanish in the limit r — > —00 . In this limit, the 
evolution of the scale factors e'^* is of Kasner type. 

(2) l^ and the physical time t all diverge in the limit r — > Tqo where 
Too is finite. In this limit, the scale factors (e"** , e^ , e'*' ) evolve as in the 
standard FRW case and (e^ , ■ ■ ■ , e'^ ) reach constant values. 

(3) Ta given in equation fl96|) provides a good approximation to t^ when 



^^ There are two special choices for the set of K^ . One is where K^ = . . . = ^^4 ^^^ 
another is the one which maximises the approximation Tq given in equation (j96|) . The 
later set may be determined by the algorithm given in Appendix D. 
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e'o ^ 1 for all / . 

(4) Any values for the ratios ^^ and ^ can be obtained, but a corre- 
sponding fine tuning of p/o = e'o is necessary. 

(5) The log periodic oscillations of l^ — l"^ , equivalently of (A\ ■ ■ ■ , A^) , 
can also be seen in the limit r — )■ Too • They can be matched to solutions of 
the type given in equation fll02p . 

To illustrate the values of Too and the ratios {^fp^, ^) one obtains, and 
to give an idea of their dependence on the initial values Iq , we tabulate these 
quantities in Table I for a few sets of initial values [Iq, K^ , U) . We have 
also tabulated the values of r^ as given by equation f l96|) . The value of gs 
follows from gg = {^^Y ^'^d, hence, is not tabulated. 





I'O' ^0' ^0' ^o) 


Ta 


T"oo 


Mil 
M4 


Ms 
Mi 


(1) 


(2,5,8,8) 


1.88 


3.21 


5.77*10-2 


2.86*10-2 


(2) 


(5,4,6,9) 


2.96 


4.16 


4.56*10-2 


1.93*10-2 


(3) 


(15,12,10,16) 


4.88 


6.61 


5.95*10-2 


1.96*10-2 


(4) 


(25,26,27,28) 


22.00 


22.59 


1.99*10-^ 


7.30*10-1° 


(5) 


(41,30,50,43) 


25.80 


28.30 


1.87*10-1° 


2.92*10-" 


(6) 


(44.5,34,49,49.5) 


34.82 


36.20 


2.59*10-1^ 


3.80*10-15 



Table I : The initial values — (/g, Iq-,Io, Iq) and the resulting values ofTa , 
Too , ^j^ , and H^ . The values in the last four columns have been rounded 
off to two decimal places. 
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In Table II, the corresponding initial values {K^,U), i = 1,2,6,7,8,9, 
are tabulated upto overall positive constants. The remaining L*s are given 
by equations fl76|) and the overall positive constants are determined hj E = 
1 and equation flTTl) . All the sets of initial values {l^^K^ ,U) are chosen 
arbitrarily with no particular pattern and are presented here to give an idea 
of the typical results. 





{K\K'^,K^,K^) oc 


{L\L'',L\L\L\V') oc 


(1) 


(4.65, 9.14, 4.57, 6.87) 


(0.60, 0.62, 0.76, 0.72,-0.94,-0.26) 


(2) 


(8.86, 8.26, 6.01, 6.62) 


-(0.08,-0.93, 0.08,-0.72, 0.54, 0.63) 


(3) 


(1.61, 2.65, 0.69, 2.1) 


-(0.2,-0.68,-0.14, 0.3, 0.08, 0.19) 


(4) 


(1.03, 1.18, 1.17, 1.27) 


(0.08, 0.58, 0.27, 0.27,-0.66,-0.66) 


(5) 


(5.24, 4.83, 4.30, 4.96) 


(0.74, 0.02, 0.24,-0.22,-0.61,-0.75) 


(6) 


(33.79, 24.23, 35.4, 32.29) 


(11.72, 9.31, 4.59,-6.46,-21.02,-21.02) 



Table II : The initial values of {K^ , U) for the data shown in Table 
I, tabulated here upto overall positive constants. These constants and the 
remaining Us are to be fixed as explained in the text. 

We find, by analysing numerous sets of initial values, that changing the 
values of {K\ U) for a given set of l^ changes the values of ^j^ and 
^ only upto about four orders of magnitude. Any bigger change requires 
changing e'o to a similar order, confirming that any values for ^j^ and ^ 
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can be obtained but only by fine tuning pjq = e'o . 

We illustrate the evolution of the universe for the data set (3) given in 
Tables I and II where many features can be seen clearly. The evolution with 
respect to r of l^ is shown in Figures 1,2 a, and 2 b. For negative values 
of r not shown in Figure 1, all l^ evolve along straight lines with no further 
crossings and their evolution is of Kasner type. Also, all l^ diverge at a finite 
value Too — 6.612 of r . The magnified plots in Figures 2 a and 2 b for 
r > 6.40 and for r > 6.55 respectively show the continually criss-crossing 
evolution of l^ which, near Too , represent the log periodic oscillations and 
are well described by equation fll02p . 
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Figure 1: The plots of l^ with respect to r . The lines continue with no 
further crossings for negative values of r not shown in the figure. All f 
diverge at t^o — 6.612 . All figures in this paper are for the data set (3) given 
in Tables I and II. 
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(a) 
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Figure 2: (a), (b) The magnified plots of l^ with respect to r for r > 6.40 
and for r > 6.55 showing the continually criss-crossing evolution of f . 
Near Too — 6.612 , these crossings are well described by equation UOS\) . 

The evolution with respect to In t of (A^, ■ ■ ■ , A^) is shown in Figure 3. 
It can be seen that (A^, ■ ■ ■ , A^) , and hence the scale factors (e'^ , ■ 
the brane directions, all stabilise to constant values as t — )■ oo . 



=A^ 



of 
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Figure 3: The plots of {X^, ■ ■ ■ , X^) with respect to In t . The lines, from 
top to bottom at the right most end, correspond to (A^, A^, A^, A^, A"^, A'^, A^) . 
(A^, ■ ■ ■ , A'^) all stabilise to constant values as t —)■ oo . 

The evolution with respect to In t of (A^, A^, \^^) and A'^ = YjI=i A* is 
shown in Figure 4. Note that the seven dimensional volume of the brane 
directions oc e^" and that it stabilises to a constant value e^" as t — )■ oo . 
We have also verified that the evolution of (A^, A^, A^'^) as t — )■ oo is same 
as that of the corresponding ones in a four dimensional radiation dominated 
FRW universe. 
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Figure 4: The plots of{X^, A^, A^°, A'^) with respect to In t . The seven dimen- 
sional volume of the brane directions oc e^" . The evolution of {X^ , X^ , X^^) 
as t —>■ oo is same as that of the corresponding ones in a four dimensional 
radiation dominated FRW universe. 



The log periodic oscillations of A'^ are illustrated in Figures 5 a and 5 b 
by magnifying the plots of {A'^—v'^) with respect to In t for In t > 20 and for 
In t > 30 . The internode seperations can be seen in these figures, and they 
match the value ~ 9.734 obtained in equation f llOSp from the asymptotic 
analysis. 
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(a) 



(b) 



Figure 5: (a), (b) The magnified plots of (A'^ — v'^) with respect to In t 
for In t > 20 and for In t > 30 showing log periodic oscillations and the 
internode seperation which is ~ 9.734. 

In all the cases we have analysed, the evolutions of {l^ , X^) are qualita- 
tively similar to the ones shown in the figures above. The details, such as the 
rise and fall of A* in the initial times or the value of Too or the stabilised val- 
ues of the brane directions, depend on the initial values but the asymptotic 
features described in the beginning of this section are all the same. Hence, 
we have presented the plots for one illustrative set of initial values only. 

VIII. Summary and Conclusions 

We summarise the main results of the paper. We assume that the early 
universe in M theory is homogeneous, anisotropic, and is dominated by A^ = 4 
mutually BPS 22'55' intersecting brane configurations which are assumed to 
be the most entropic ones. Also, the ten dimensional space is assumed to be 
toroidal. We further assumed that the brane antibrane annihilation effects 
are negligible during the evolution of the universe atleast until the brane 
directions are stabilised resulting in an effective 3 + 1 dimensional universe. 

We then present a thorough analysis of the evolution of such an universe. 
We obtain general relations among the components of the energy momentum 
tensor T45 using U duality symmetries of M theory and show that these 
relations alone imply, under a technical assumption, that the A^ = 4 mutually 
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BPS 22'55' intersecting brane configurations with identical numbers of branes 
and antibranes will asymptotically lead to an effective (3 + 1) dimensional 
expanding universe. 

To obtain further details of the evolution, we make further assumptions 
about Tab ■ We then analyse the evolution equations in D dimensions in 
general, and then specialise to the eleven dimensional case of interest here. 
Since explicit solutions are not available, we apply the general analysis and 
describe the qualitative features of the evolution of the A^ = 4 brane con- 
figuration : In the asymptotic limit, three spatial directions expand as in 
the standard FRW universe and the remaining seven spatial directions reach 
constant, stabilised values. These values depend on the initial conditions and 
can be obtained numerically. Also, any stabilised values may be obtained but 
it requires a fine tuning of the initial brane densities. 

We also present a physical description of the mechanism of stabilisation 
of the seven brane directions. The stabilisation is due, in essence, to the 
relations among the components of Tab which follow from U duality sym- 
metries, and to each of the brane directions in the A^ = 4 configuration 
being wrapped by, and being transverse to, just the right number and kind 
of branes. This mechanism is very different from the ones proposed in string 
theory or in brane gas models. 

In the asymptotic limit, from the perspective of four dimensional space- 
time, we obtain an effective four dimensional Newton's constant 6*4 which is 
now time varying. Its precise time dependence will follow from explicit solu- 
tions of the eleven dimensional evolution equations. We find that, in the case 
of A^ = 4 brane configuration, G4 has characteristic log periodic oscillations. 
The oscillation 'period' depends only on the configuration parameters. 

Using numerical analysis, we have confirmed the qualitative features men- 
tioned above. 

We now make a few comments on the assumptions made in this paper. 
Note that the assumptions mentioned above in the first paragraph of this 
section pull a rug over many important dynamical questions that must be 
answered in a final analysis. Some of these questions, lj in the context of M 
theory, are: 

* Starting from the highly energetic and highly interacting M theory 
excitations, which are expected to describe the high temperature state of the 
universe, how does a eleven dimensional spacetime emerge? 



^•^Many of the questions listed below have been raised by the referee also. 
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* What determines the topology of the ten dimensional space? Here, we 
assumed it to be toroidal. How does the universe evolve if its spatial topology 
is not toroidal? 

* From what stage onwards, is the eleven dimensional 'low energy' effec- 
tive action a good description of further evolution? 

* What are the relevant 'low energy' configurations of M theory? Here, 
based on the black hole studies, we have assumed that the A^ = 4 mutually 
BPS 22'55' intersecting brane configurations are the most entropic ones and, 
hence, that they are the dominant configurations in the early universe studied 
here. 

This raises further questions: Are the 22'55' , and not some other mu- 
tually BPS A^ > 4 or some other non BPS, configurations really the most 
entropic and the dominant ones? Even assuming that mutually BPS A^ = 4 
is the answer, are there other A^ = 4 configurations beside the 22'55' ones 
and, if so, how do they affect the evolution described here? What are the 
effects of the subdominant configurations? In particular, will the effects of 
other brane configurations mentioned above undo the stabilisation of seven 
directions presented here? 

Note that unless these questions are answered and, furthermore, it is 
shown that other brane configurations mentioned above do not undo the 
stabilisation presented here, our assumption that the evolution of the uni- 
verse is dictated by the 22'55' configuration amounts to a fine tuning: The 
22'55' configuration assumed here, where the sets of 2 branes and 5 branes 
wrap the directions {x^, • ■ ■ , x'^) homogeneously everywhere in the mutually 
transverse three dimensional space, may not arise generically. Also, the im- 
plicitly required absence of other brane configurations is not natural in the 
context of early universe. Then the problem of the emergence of an effec- 
tive 3 + 1 dimensional universe, a solution for which is presented here, gets 
shifted to answering how the required, finely tuned, initial conditions may 
arise naturally from M theory. 

* What is the time scale of brane antibrane annihilations in the 22'55' 
configuration studied here? Is it long enough for the brane directions to be 
stabilised as described in this paper? Here, based on the black hole studies, 
we have assumed it to be long enough. 

* A related question, but applicable after stabilisation of brane directions, 
is the following: If all the branes and antibranes will eventually decay, as 
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seems natural, then what are the decay products? How can one obtain the 
known constituents of our present universe? 

Although one of us have presented a principle in [8] that may be of help, 
the fact is that we do not know even where to begin in answering these ques- 
tions quantitatively, much less know the answers. Nevertheless we present 
the above list of questions, unlikely to be complete, in order to emphasise 
the further work required to understand how our known 3 + 1 dimensional 
universe may emerge from M theory. 

In the present work, with many attendant assumptions, we considered the 
22'55' configurations and explained a mechanism by which seven directions 
stabilise and an effective 3 + 1 dimensional universe results. Clearly, it is 
important to answer the questions listed above and thereby determine the 
relevance of this mechanism. 

Within the present framework, there are many other issues that may be 
studied further. We conclude by mentioning a sample of them. We have 
shown here that a large stabilised seven dimensional volume can be obtained 
but it requires a corresponding fine tuning of initial brane densities. This is 
within the context of our ansatzes for T^b and the equations of state. It 
will be of interest to prove or disprove the necessity of such a fine tuning in 
more general contexts. 

The A^ = 4 intersecting brane configuration studied here is the entropi- 
cally favourable one and, as proposed in [S] , may be thought of as emerging 
from the high temperature phase of M theory in the early universe. Such an 
emergence suggests that there may be novel solutions to the horizon problem 
and to the primordial density fluctuations, perhaps similar to those explored 
recently in the Hagedorn phase of string theory by Nayeri et al [52]. Note 
that this involves answering many of the questions listed above. 

It may be of interest to study further the consequences of time varying 
Newton's constant which appears here, in particular possible imprints of its 
asymptotic log periodic oscillations. 

In the case of a class of black holes, the brane configurations describe 
well their entropy and Hawking radiation. In the present description of a 
four dimensional early universe in terms of A^ = 4 intersecting branes, it 
is not clear which quantities to calculate which, analogously to entropy or 
Hawking radiation in the black hole case, may provide further validation. It 
is important to study this further. 

Acknowledgements: We thank Sanatan Digal for help in numerical 
analysis, and are grateful to the referee for her/his helpful comments. 
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Appendix A : U duality relations in black hole case 

Consider black holes in m + 2 dimensional spacetime described by mu- 
tually BPS intersecting brane configurations in M theory. The brane action 
Sbr in equation (|3]) is the standard one for higher form gauge fields. The 
corresponding black hole solutions and their properties are well known, so 
here we only highlight the points related to U duality symmetries. Also, for 
illustration, we consider only 2 branes and 5 branes. 

As mentioned in section II, the method of U duality symmetries applies 
here also and leads to the same relations between A* . They are best seen 
in the extremal case. (The non extremal case requires further analysis and 
is more involved.) The eleven dimensional line element ds for the extremal 
brane configurations are of the form 

ds^ = -e^^°dt^ + J2e'^\dx'y (105) 

i 

where (A°, A*) depend on r, the radial coordinate of the m + 1 dimensional 
transverse space. For 2 branes and 5 branes, the A*s may be written as 

X' = X' = -'^ , X^ = ... = X^^ = l (106) 



Ai° = 


h 
^ 6 


Ai° = 


2h 
6 



V = --- = A^ = -^ , A^ = ... = A^° = ^ (107) 

where e^ = H = 1 + -^^ is the corresponding harmonic function and Q is 
the charge. See, for example, [53] for more details. Clearly, the U duality 
relations in equations flTTj) . f lT2|) . and f lT3|) are valid here also. 

For the extremal 22'55' configuration (12, 34, 13567, 24567) , the trans- 
verse space is three dimensional and the A*s may be written as 

A' = - {-2hi + h2-h + 2/^4) 
6 

A^ = - (-2/ii + h2 + 2h3 - hi) 
6 

A^ = - (hi - 2h2 -h + 2hi) 
6 

A^ = ^ {h - 2h2 + 2h - K) 
6 
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A^ = A^ = A^° = lihi + h2 + 2h + 2K) (108) 

o 

where e'^' = H^ = 1 + — are the corresponding harmonic functions and 
Q^s are the charges. Clearly, the U duality relations in equations f lT^ and 
f lT5]) are valid here also. Furthermore, if 2 and 2' branes are identical then 
hi = /i2 and we get A^ = A^ , and similarly other relations when different 
sets of branes are identical. 

We further illustrate the U duality methods by interpreting a U dual- 
ity relation J2i ^iX = as implying a relation among the components of 
the energy momentum tensor T^^ . The relations thus obtained are indeed 
obeyed by the components of Tab calculated explicitly using the correspond- 
ing higher form gauge field action S^t ■ 

For this purpose, let the spacetime coordinates be x^ = (r, a;°) where 
x" = (x°, x\ 6"") with x^ = t , i = 1,- ■ ■ ,q , a = l,---,m, and q + m = 9 . 
The x* directions may be taken to be toroidal, some or all of which are 
wrapped by branes, and 9"' are coordinates for an m dimensional space of 
constant curvature given by e = ±1 or . The metric and brane fields depend 
only on r coordinate. We write the line element ds, in an obvious notation, 
as 

ds'' = -e-'^odf + ^ e^^'idx^ + e^Hr' + e^^^d^l^^ . (109) 

i 

The independent non vanishing components of T\ are given by T''^ = / 
and T^Q, = p^ where a = {0,i,a) . These components can be calculated 
explicitly using the action S^r ■ For example, for an electric p-brane along 
{x^,- ■ ■ , x^) directions, they are given by 

PO = Pll = -Pi. = -Pa = f = ^^^^ F01...p. F°l-^'- (110) 

where p\\ = pi ior i = 1,- ■ ■ ,p , p± = pi ior i = p + 1, ■ ■ ■ , q , and note 
that / is negative. For mutually BPS A^ intersecting brane configurations, 
it turns out [32] - [51] that the respective energy momentum tensors T"^ 
and T^Bii) obey equations analogous to those given in fl2T|) - f l23|) . 

Equations of motion may now be obtained from equations (jl]) and fl2T]) . 
After some manipulations, they may be written in a form similar to those 
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given in d?]) - (jH]) as follows: 

A'-E(^r)' = 2/ + em(m-l)e-2- (111) 

a 

+ e (m - 1)6-^" 5"" (112) 

fr + fAr-T.Po^K = (113) 

a 

where A = J2a ^"^ = -^°+I]i A*+mcr and the subscripts r denote r-derivatives. 
See [50j particulary, whose set up and the equations of motions are closest 
to the present ones. 

Consider now the case of 2 branes or 5 branes. We assume that pa = p± 
which is natural since 6°- directions are transverse to the branes. Applying 
the U duality relations in equation f lT3|) then implies, for both 2 branes and 
5 branes, the relation 

P\\=Po+P± + f (114) 

among the components of their energy momentum tensors. Note that it 
is also natural to take po = P\\ since x^ = t is one of the worldvolume 
coordinates and may naturally be taken to be on the same footing as the 
other ones (x^, ■ ■ ■ , x-^) . Equation (11141) then implies that p± = —f . The 
relation between py and / is to be specified by an equation of state which, 
in the black hole case, is that given in equation (1 11 01) . 

For now, however, we take Pq and py to be different. Keeping in mind 
that / is negative, we assume the equations of state to be of the form pai = 
— (1 — n^) // where a = {0,i,a) , I = 1,- ■ ■ ,N , and n^ are constants. Then 
for 2 branes and 5 branes, we have 
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2 : Uo, = {uo, u\\, ui\, uj_, uj_, u_i_, Uj_, u_i_, m_l, Uj_) 

5 : Ua = {uo, uii, u\\, m, u\\, M||, u±, u±, u±, u±) (115) 

where the / superscripts have been omitted here and My = uq + u± which 
follows from equation (I114p . Note that -u^ = and uq = uy = 2 in the 
black hole case given in equation (IllOp . 
Let 

a 

G„/3 = 1 - 6^^ , ^" = E G"-^ < ^i ■ (116) 
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Then, after a straightforward algebra, one obtains 

ll = -Y^ giJ e" + u^ e m{m - 1) e^^^-'^^ , (117) 

J 

which are similar to equations ( 139|) . The remaining equations are not needed 
and, hence, not given here. Using equations flllSp and f lll6p . it is now 
straightforward to calculate Q^'^ for A^ intersecting brane configurations. 
It turns out because of the BPS intersection rules that Q^'^ may be written 
as 

G^^ = 2uo {u± - uo5^^) . (118) 

The corresponding Qu is given by 



1 . u^ 
2uq Nu± — uq 



Qij = —A .r.. .. -M • (119) 



Now take po = P\\ ■ Then equation (I114p gives p± + / = . In terms of Ua , 
we now have uq = u\\ and m_l = . Clearly, then Q^'^ oc 6^'^ and equations 
(11171) can be solved for /^(r). See [50] for such solutions, with uq = 2 as 
follows from equation (IllOp . and their analysis. 

Tracing through the steps in the above derivation, it can also be seen that 
for the homogeneous early universe case, (r, /) here get replaced by {t, —p) , 
and (t,pQ,UQ) here get replaced by {r,p±,u± = u) . Then, equations (I114p 
and (IllSp become equations ( l20l) with z = —1 and ( 172|) respectively. 



Appendix B : To show E >Q 

Let 1 = (1, 1, ■ ■ ■ , 1) and v = {vi, V2,- ■ ■ , Vn) be the standard n - 
component vectors with the standard vector product. Let 9n be the angle 
between them. Then 1 ■ 1 = n , v ■ v = Z^^'y^ ' (1 ' ^Y — ij^a'^^a)'^ = 
n cos^On X^a ^a ' ^^'1 ^^ have the Schwarz inequality in the form 



n n 

n 

a=l a=l 



T.<-ii: ^af = nal>0 (120) 



where a^ = sin^On Yl^=i ^1 ■ The equality is valid, i.e. (T„ = , if and only 
if sin 6n = , equivalently Vi = ■ ■ ■ = Vn ■ 
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We now show the following: 

Let C^ and Gij be given by equation (^ - Ifui and D satisfy the relations 
J2i UiU = and Y.ij G'^UiUj > then 2E = - Y.ij GijDU >0 . E vanishes 
if and only if U all vanish. 

Proof : It is clear that E vanishes if U all vanish. Now, let 1 = 
(1, 1, ■ ■ ■ , 1) , M = {ui, ■■■,ud^i) , and 9 be the angle between them. 
Then {EiUi^ = iD-l) cos^O Y.iu\ . Hence, Y^aG'^iU^ = ^ (Ei^i)' - 
Y.i'^'i > implies that 

1- (D- 1) sin^e > . (121) 

The vector L = (L^, ■ ■ ■ , L^~^) is perpendicular to u since J^iUiL^ = . 
Let L = L± + L|| where L± is perpendicular to the plane defined by 1 and 
u , and L|| lies in it. Then I]j(-^*)^ = -^1 + -^y where L^ = L± ■ L± and 
Lm = L|| ■ L|| . Since L and H are perpendicular and Ly lies in the plane 
defined by 1 and u , it follows that Ly is perpendicular to u, and that the 
angle between the vectors 1 and Ly is | ± ^ . We then have 



2E = -Y.G,,VU=Y.{Vf-{Y.Lr 

ij i i 

= Ll + Ll-{D-1) L\ sin^e > 

where the inequality follows from equation f ll2ip . The equality holds, and 
hence E vanishes, only when L^ = L| = , z.e. only when L* all vanish. 
This completes the proof. 



Appendix C : Signs and non vanishing of (A,-, l!^) 



Here, we show that the inequality in equation (ISO]) implies that none of 
(At-, Z:^) may vanish, and that they must all have same sign. 

Setting xi = li , equation ([80]) becomes X = 12u'^ {E + J2ie^') > 
where the polynomial X = {X1+X2 + X3 + ^4)^ — 3(xf + xl + xl + xl) . Now, 
if any of the xj vanishes then X < , see the Scwarz inequality given in 
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equation (11201) . Hence, none of the xj may vanish. Rewrite X as 

X = \^{Xi + X2 + X^y -3{xl + xl + xDj -2xl + 2Xi {Xi + X2 + X3) 

= [{xi + X2f - 2{xl + xl)] + {(X3 + Xif - 2(x^ + xD] 

-{xf + xl + xl + xl) + 2(xi + X2) (X3 + X4) 

and note that {■■■}< for each curly bracket, see equation f ll20p . Hence, 
the necessary conditions for X > are 

X4 {xi + X2 + 2:3) > , (xi + X2) (2:3 + 0:4) > . 

Let one of the x/ , e.g. X4 , be negative and the other three positive. This 
violates the first inequality above and, hence, is not possible. Let two of the 
xi , e.g. X3 and X4 , be negative and the other two positive. This violates 
the second inequality above and, hence, is not possible. Similarly, three of 
the xj being negative and one positive is also not possible. Thus, the only 
possibility is that all xj have same sign. Thus we have that none of the l!;. 
may vanish, and that they must all have same sign. 

With l^ denoted as xi , equation f l79|l for A^ becomes 

6u At- = 2xi + 2x2 + X3 + X4 + 6m L . 

Note that u > . If L = then it follows that A,- does not vanish and has 
the same sign as x/ . Consider now the case where L 7^ . Using equation 
f ll2UI) to eliminate Y^jxj in the polynomial X, we obtain 

X = - (xi + X2 + X3 + X4)^ - 3al = 12^2 {E + J2 e'') • 

Using the inequality 2E > 3(L)^ , see equation (178ll . it follows that (xi+X2 + 
X3 + X4)^ > 72m^(L)^ . Combined with the earlier result on l!^ , this inequality 
implies that (xi + X2 + X3 + X4 + QuL) , and hence A,- given above, may not 
vanish and must have the same sign as xj = ll , irrespective of whether L is 
positive or negative. This completes the proof. 



Appendix D : Set of K^ which maximises r^ 
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With no loss of generality, let < — /q < • ■ ■ < —Iq . The corresponding 

set of K^ which satisfies equation fl7ip . with E = 1 , and which maximises 

ii 
Ta = min{Ti} , where tj = —j^ , may be obtained by the following al- 
gorithm. The required analysis is straightforward but a little tedious and, 
hence, is omitted. 

• Let K^ = — /q K . It will turn out that Tq = ri = -^ . 

• Choose K^ = -llK . Then ra = n . 

• If -ll -ll < -ll then choose K^ = K^ = -{l^ + l^) K . Then 

T4 > 7-3 > r2 = Ti . 

• If -ll - /q > -^0 then choose K^ = -If, K . Then Ts = r2 = ri . 

• If -ll -ll > -ll and if -ll - ll - tl < -21^ then choose K^ = 
-\ ill + ll + ll)K . Then T4 > rg = T2 = n . 

• If -ll -ll> -ll and if -ll -ll-ll> -211 then choose K^ = -It K . 
Then r4 = ■ • ■ = Ti . 

• K^ are all thus determined in terms of K . Equation fl7T|) . with £" = 1 , 
will now determine K . 
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